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Abstract 

The solar rotation profile is well constrained down to about 0.25 Rq thanks to 
the study of acoustic modes. Since the radius of the inner turning point of a resonant 
acoustic mode is inversely proportional to the ratio of its frequency to its degree, 
only the low-degree p modes reach the core. The higher the order of these modes, the 
deeper they penetrate into the Sun and thus they carry more diagnostic information 
on the inner regions. Unfortunately, the estimates of frequency splittings at high 
frequency from Sun-as-a-star measurements have higher observational errors due to 
mode blending, resulting in weaker constraints on the rotation profile in the inner 
core. Therefore inversions for the solar internal rotation use only modes below 2.4 
mHz for ^ < 3. In the work presented here, we used an 11.5 year-long time series to 
compute the rotational frequency splittings for modes < 3 using velocities measured 
with the GOLF instrument. We carried out a theoretical study of the influence of 
the low-degree modes in the region 2 to 3.5 mHz on the inferred rotation profile as 
a function of their error bars. 

Keywords: Helioseismology, Observations, Inverse Modeling; Interior, Radiative 
zone. Core; Rotation 
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1. Introduction 

Our knowledge of the solar rotation profile has been derived from the study of the 
resonant acoustic modes, which are trapped in the solar interior. Since the solar ro- 
tation lifts the azimuthal degeneracy of these resonant modes, their eigenfrequencies 
{Vnim) split into their m-components: where £ is the angular degree, n the radial 
order, and m the azimuthal order. This separation Au^i^ — usually called rotational 
splitting (or just splitting) — depends on the rotation rate in the region sampled by 
the mode. Using inversion techniques, the rotation rate at different locations inside 
the Sun can be inferred from a suitable linear combination of the measured rotational 
splittings 

Today, the rotation rate inside the Sun is rather well known above 0.4 Rq 
(Thompson et al, 1996; Schou et al, 1998; Howe et al, 2000; Antia and Basu, 
2000). The convective zone is characterized by a differential rotation extending from 
the surface down to the tachocline, located around 0.7 Rq. Below the tachocline, 
inside the radiative region, the Sun appears to rotate as a rigid body with a nearly 
constant rate of f» 433 uHz down to the solar core, i.e. « 0.25 Rq. The rotation rate 
inside the core derived from p modes is still uncertain (Jimenez et al., 1994; Elsworth 
et al, 1995; Chaplin et al., 2001; Chaplin et al, 2004). Recent measurements of the 
asymptotic properties of the dipole g modes and their comparison with solar models 
favors a faster rotation rate inside the solar core (Garcia et al., 2007). 

For p modes of a given degree i, the radius at the inner turning point is a 
decreasing function of frequency, given by: 

n = Ct L/ {2'KVnl) 

where L = I -\- 1/2, Vnf, is the frequency of the mode, and ct = c(rt) the sound- 
speed at the radius rt (see for example. Lopes and Turck-Chieze, 1994). Thus the 
modes with increasing frequencies - higher radial order n - penetrate deeper inside 
the Sun. Unfortunately, when fitting Sun-as-a-star observations the uncertainties 
on the rotational splittings that are the most sensitive to low-degree p modes are 
very large. Indeed, as the modes lifetimes decrease with frequency their line widths 
increase. Therefore, for frequencies above ~ 2.3 niHz, there is a substantial blending 
between the visible m-components of the p modes. This blending makes it difficult 
to extract precisely the rotational splitting. At higher frequencies (« 3.9 mHz) even 
the successive pairs of £ = 0,2 and £ = 1,3 modes blend together. As a result, 
with today's fitting methods, it is still not possible to obtain values of the rotational 
splittings with an accuracy good enough to be useful in any rotation inversion. 

By contrast, at low frequency - below n = 16, or about 2.4 mHz - the lifetime 
of the modes increases and thus their line width is very small. This allows us to 
measure their rotational splittings with very high precision. However, these modes 
have inner turning points at shallower depths than the high-frequency modes (above 
0.08 and 0.12 Rq for the 1=1 and 2 modes respectively). Therefore, even though 
these modes do not carry any information below » 0.1 Rq, they help improve our 
knowledge of the inner rotation rate because their inclusion contributes to an increase 
of the precision of the inversions since they have smaller error bars (Eff-Darwich, 
Korzennik, and Jimenez-Reyes, 2002; Eff-Darwich et al., 2007). 

For all of these reasons, rotation inversion methodologies usually limit the input 
data set of low-degree p modes to low frequency modes. For example, Couvidat et al. 
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(2003) used a limited number of splittings, corresponding to modes with £ < 3 and 
frequencies below 2.4 mHz (n = 15) resulting from fitting GOLF^ and MDI^ 2,243 
day-long velocity time series, to infer the solar rotation profile. He concluded that the 
uncertainties in the rotation rate below 0.3Rq were still quite large. Therefore, to 
obtain a better and more reliable rotation profile in the inner core, we need, on one 
hand, to include low-frequency acoustic, mixed, and gravity modes in the inversions 
(see discussions in Provost, Berthomieu, and Morel, 2000 and Mathur et al, 2007). 
On the other hand, we also need to measure more accurately low-degree high-order 
p modes to push further the frequency limit of the modes used in the inversions. 

We present, in Section 2, rotational splittings of / < 3 modes computed using 11.5 
years of GOLF data, paying special attention that no bias is introduced by the solar 
activity or by the length of the fitting window. Thanks to the length of these new 
time series, the splitting error bars are improved by y/T when compared to previous 
and shorter analysis using the same instrument. In Section 3, we study the sensitivity 
of the rotation rate below 0.25 Rq to the available p modes and their error bars. 
We then discuss, in Section 4, the influence of these modes on the inversions, based 
on the resulting resolution kernels and we flnish by discussing the inverted rotation 
profile using the newly computed GOLF splittings. 



2. Data Analysis 

We analyzed a 4 182 day-long time series of GOLF observations that spans the 11 
April 1996 to 22 September 2007 epoch (w 11.5 years). These observations were cali- 
brated following the methods described in Garcia et al. (2005). The duty cycle of that 
time series is 94.5%. The GOLF instrument is a resonant scattering spectrophotome- 
ter onboard the SOHO^ spacecraft. A standard Fast Fourier Transform algorithm 
was used to compute the power spectral density. This periodogram estimator was 
used because our primary interest is the study of modes in the medium frequency 
range. We have therefore decided to privilege the simplest estimator, avoiding the 
use of multi-tapers or zero-padding estimators that are better suited to look for weak 
and narrow peaks at low frequency. 

To obtain the mode parameters, each pair of modes {£ = 0,2 and £ = 1,3) 
was fitted to a set of asymmetrical Lorentzian, as defined in Nigam and Kosovichev 
(1998), using a maximum- likelihood method since the power-spectrum estimator 
follows a distribution with two degrees of freedom. For additional details on this 
method see, for example, Toutain and Appourchaux (1994). For each mode, we fitted 
the following set of parameters: the amplitude, the line width, the central frequency, 
the peak asymmetry, the rotational splitting (if £ > 0), and the background noise. 

It has been shown that the amplitude ratios of the m components inside a multi- 
plet could produce a bias in the splitting determination if they are not correctly set 
(Chaplin et al, 2006). We have therefore fixed them to the averaged values obtained 
directly from the GOLF data set itself (for details see Henney, 1999). Finally, to 
reduce the parameter space, we used the same asymmetry and line width for both 



Global Oscillations at Low Frequency (Gabriel et al., 1995) 
■^Michelson Doppler Imager (Scherrer et al., 1995) 

^ solar and Heliospheric Observatory (Domingo, Fleck, and Poland, 1995) 
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modes in any given fitting window. Indeed Chaplin et al. (2006) demonstrated that 
using the same width for the modes fitted simuhaneously in a given fitting window 
reduces the bias in the splitting determination at high frequencies. Finally, the un- 
certainties on the fitted parameters were derived from the square root of the diagonal 
elements of the inverted Hessian matrix. 




•= 400 
^ 350 



1500 2000 2500 3000 3500 1500 2000 2500 3000 3500 
Frequency {fitiz) Frequency (fiHz) 

Figure 1. Synodic splittings for the £ = 1,2, and 3 modes fitted using increasing fitting 
window widths: from 30 to 65 /iRz. The vertical dotted lines mark the upper limit of the 
modes whose splittings are used in this work. The error bar in the bottom-left corner of each 
plot corresponds to ± 20 nHz and is plotted as reference. 



It has been shown, using artificial simulations, that the determination of the 
rotational splitting could be biased by the width of the fitting window (Chaplin 
et at, 2006). This is caused by the leakage of the neighboring modes. For example the 
splittings of the £ = 2 modes could be noticeably modified (increased) at frequencies 
above 2.5 mHz due to the leakage of the £ = 4 modes. We have studied this effect 
using the GOLF data by selecting 4 different windows for the fit: 30, 45, 55, and 
65 /iHz. The results are plotted in Figure 1. As expected, the values fitted for the 
£ = 2 splittings increase at high frequencies when wider fitting windows are used. 
This behavior is consistent with the results obtained by the solarFLAG group (see 
for comparison Figure 9 of Chaplin et al, 2006). However, the smallest window (30 
/xHz), also biases the fittings at high frequency because fitting such narrow spectral 
range prevents the fitting code from properly constraining the background level. 
Therefore, we have adopted the 45 /iHz window as the best compromise between 
these two effects. In the case of the £ = I modes the splittings computed using 
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Figure 2. (a) Fitted splittings for the mode {£ = 2,n = 17) as a function of the solar activity 
level measured by the solar radio fiux at 10.7cm (Isfu = 10~^^Wm~^Hz~^). The asterisk 
writh thick error bar indicates the result of fitting the full time series, (b) Sensitivity of the 
splittings to the activity for all of the fitted modes (£ = 1,2 and 3) . 



the smaller windows are roughly constant up to 3.4 mHz and then increase. For 
higher frequencies the sphttings increase whatever window width is used and it 
is not possible to distinguish between a bias induced by the blending of the two 
adjacent m-components or a real increase. Therefore we have hmited our study to 
the splittings of modes below 3.4 mHz (this limit is indicated with a vertical dotted 
line in Figure 1) where all of the splittings are roughly constant within the error 
bars. Another approach to limit the leakage of the neighbor peaks could be to fit the 
whole spectrum at once with a simple model of the background noise (Roca Cortes 
et ai, 1998). 

It is also important to notice that some splittings (for example the ^ = 1 n = 17 
mode) are slightly out of the general trend of the rest of the splittings (less than 
2 (t) as a consequence of the stochastic nature of the excitation. This deviation is 
in quantitative agreement with the results obtained from artificial simulations (see 
Chaplin et ai, 2006). 

We have also verified that there is no noticeable influence of the solar activity 
on the extracted splittings. To do this, we followed an approach similar to the one 
described in Jimenez-Reyes et al. (2001): we divided the time series in 12 independent 
350 day-long segments'* and fitted them simultaneously. A linear dependence on a 
solar activity index was added to each parameter defining the spectral profiles, with 
two exceptions: first, we fixed the asymmetries to the values found by analyzing 
the full time series; second, we included a quadratic dependence on activity for the 
frequencies of the modes. We used the integrated 10.7 cm radio fiux measurements 
from the National Geophysical Data Center"', averaged over the duration of each 
sub-series as indicator of solar activity. The 10.7 cm radio flux has been shown to 
be the indicator that best correlates with the measured changes resulting from solar 
activity. 

Figure 2a illustrates the variation of the splitting with solar activity for one mode 
{I = 2, n = 17), while Figure 2b shows the sensitivity of the splittings to the radio 



*one 350 day-long segment was not fitted since its duty cycle is very small because of the loss 
of contact with the SOHO spacecraft during that epoch 

■''http: / /www. ngdc.noaa.gov/ngdc. htm I 
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flux for aU of the fitted modes (i.e. the slope in Figure 2a). This plot shows that there 
is no clear systematic dependence like the one seen for the frequency (Jimenez, Roca 
Cortes, and Jimenez-Reyes, 2002; Gelly et at, 2002; Garcia et al, 2004). Moreover 
variations during the solar cycle are generally marginally significant. The splittings 
obtained by fitting the whole time series are fully consistent with those obtained 
when using this activity-dependent method and correspond to the Sun at its mean 
activity level, as expected. The resulting rotational splittings for modes £ < 3 are 
listed in Table 1. 

Before inverting these rotational frequency splittings we completed the set of 
splittings of low-degree low-order p modes down to 1 mHz with those extracted from 
the analysis of combined GOLF and MDI time series (Garcia et at, 2004). Then 
we added high-degree modes using splittings of modes between £ = 4 and £ = 25 
from the analysis of the 2 088 day-long time series of MDI observations fitted by 
Korzennik (2005). 



3. Sensitivity of the Splittings to the Rotation Rate below 0.25 Rq 

We have shown in the introduction that the acoustic modes of higher radial order 
penetrate deeper in the solar interior and that they are potentially of great interest 
to better constrain the rotation inside the solar core. This is illustrated in Figure 3 
where the modes listed in Table 1 are plotted as a function of the radius of their 
inner turning point. 
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Figure 3. Fitted modes listed in Table 1 as a function of the inner turning point radii. The 
error bars arc the sphtting error bars in nHz magnified a factor 10'*. 



The vertical error bars correspond to the 1 a error bars of the splittings. These 
uncertainties are smaller for the £ = 2 and 3 than for the £ = 1 modes because 
the visible m components of the latter are closer in frequency and therefore blend 
together at lower frequencies than for higher-degree modes. 
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Table 1. Central frequencies {uo), synodic sectoral 
splittings [Au^g^), and their respective 1 a error 
bars {(To and (7„ern) computed using 4182 days of 
GOLF velocity time series. 



e 


n 














(/.Hz) 


(nHz) 


(nHz) 


(nHz) 


1 


7 


1185.5893 


4.0 


402.2 


3.2 




8 


1329.6368 


3.3 


405.2 


2.9 




9 


1472.8475 


5.0 


400.0 


4.2 




10 


1612.7268 


9.9 


407.6 


8.0 




11 


1749.2887 


10.5 


388.1 


9.2 




12 


1885.0853 


12.2 


400.4 


10.6 




13 


2020.8228 


14.1 


416.5 


12.5 




14 


2156.8158 


18.5 


397.3 


18.4 




15 


2292.0340 


19.0 


404.7 


21.6 




16 


2425.6381 


20.0 


396.0 


26.3 




17 


2559.2441 


19.3 


345.2 


30.7 




18 


2693.4385 


18.5 


381.5 


27.0 




19 


2828.2585 


18.0 


387.1 


26.0 




20 


2963.4225 


17.7 


403.8 


24.9 




21 


3098.2937 


19.0 


405.8 


29.3 




22 


3233.2855 


22.2 


419.2 


41.6 




23 


3368.6923 


26.9 


385.1 


75.6 


2 


8 


1394.6851 


14.2 


401.8 


5.9 




9 


1535.8642 


7.0 


403.1 


4.1 




10 


1674.5434 


12.3 


406.5 


6.9 




11 


1810.3293 


13.3 


403.0 


6.9 




12 


1945.8173 


15.9 


391.6 


8.4 




13 


2082.1175 


18.5 


379.2 


8.6 




14 


2217.6968 


22.2 


387.4 


11.3 




15 


2352.2626 


23.7 


393.4 


13.5 




16 


2485.9310 


22.8 


414.7 


12.9 




17 


2619.7165 


20.6 


417.6 


11.7 




18 


2754.5761 


20.4 


413.5 


12.2 




19 


2889.6857 


20.5 


402.7 


12.7 




20 


3024.8519 


20.3 


443.9 


11.7 




21 


3159.9752 


22.3 


437.2 


13.6 




22 


3295.2381 


28.1 


392.3 


21.4 




23 


3430.9280 


38.6 


447.5 


33.4 


3 


9 


1591.4891 


19.5 


382.3 


7.1 




10 


1729.1460 


47.7 


417.5 


16.3 




11 


1865.3041 


36.9 


411.6 


12.9 




12 


2001.2566 


28.6 


404.9 


9.4 




13 


2137.7634 


38.2 


383.3 


11.4 




14 


2273.4932 


32.7 


407.4 


10.6 




15 


2407.6669 


36.2 


403.8 


11.5 




16 


2541.7150 


32.7 


414.1 


10.4 




17 


2676.2540 


27.1 


403.7 


8.9 




18 


2811.4796 


24.1 


397.5 


8.3 




19 


2947.1109 


23.1 


411.6 


7.8 




20 


3082.4893 


27.2 


409.9 


9.3 




21 


3217.9177 


31.3 


415.0 


11.4 




22 


3353.6765 


46.1 


388.6 


19.0 




23 


3489.6852 


66.4 


425.2 


30.0 
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The amount of information on the rotation rate inside the core (r < 0.25 i?©) 
present in these modes depends on tire precision of the measured splittings. The 
sensitivity of tire splittings to the rotation rate below 0.25 Rq is plotted in Figure 4a. 
This sensitivity is defined as the ratio of the contribution of the rotation rate below 
0.25 Rq to the splitting, ( Aungm) and the imcertainty ((t„^„j) in its determination 
(see equation 1 for sectoral modes, £ = m). 



Sensitivity 



r 

Jo 



0.25R(: 



(1) 



For a given mode, if the sensitivity function is greater than one, this particular 
mode provides useful information on the core rotation. Up to w 3.4 mHz, all of 
the modes I < 3 are potentially interesting for the inversions. The situation has 
evolved drastically since 2003 when the previous analysis of 2 243 day-long GOLF 
and MDI time series showed that below 0.2 Rq only those modes below 2.2 mHz had 
a sensitivity greater than one (see left panel of Figure 1 in Couvidat et al, 2003). 
In fact, when we compute the sensitivity function below the latter radius using the 
new computed GOLF splittings we obtain the same dependence with the frequency 
below 0.25 Rq but rescaled to lower sensitivity. For the £ = 2 and 3 modes the 
sensitivity reduction is about a factor of two. In the case of the £ = I modes the 
factor is smaller and the sensitivity curve crosses unity at 3.1 instead of 3.3 mHz. 

Figure 4b shows the resulting splittings normalized by their error bars. The pre- 
cision clearly increases towards the low-frequency range for the £ = 1 modes while 
for the £ = 2 and 3 modes the curves are flatter indicating that the error bars are 
nearly constant all along the analyzed frequency range. 
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Figure 4. (a) Sensitivity of the splittings to the rotation rate below 0.25 -Rq as defined in the 
text. The horizontal dotted line is the lower limit below which the splittings would not carry 
any useful information on the solar core rotation, (b) Splittings normalized by their error bars. 



4. Discussion 



To study the effect of the / < 3 modes at medium-range frequencies on the inferred 
rotation rate we carry out some 2D inversions. 
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The rotational splittings {Avnim) are the integral of the product of a sensitivity 
function - or a kernel - Knim{i^7(^)7 ^ known function based on a solar model, with 
the rotation rate {Q{r,9)) over the radius (r) and the co-latitude [9) (Hansen, Cox, 
and van Horn, 1977): 

At-nfm = / / Knirnir, 0)fl{r, e)drd9 + Enlm (2) 

27rJo Jo 

where e„^,„ is the effective error on the measured value Ai'nim- We assume that the 
errors (e„^m) follow a normal distribution with a standard deviation cr„fmi which is 
estimated by the fitting procedure (see Section 2), as listed in Table 1. This set of 
equations defines a classical inverse problem for the solar rotation. The inversion of 
this set of AI integral ecjuations - one for each measured Avnim ^ allows us to infer 
the rotation rate profile as a function of radius and latitude from a set of observed 
splittings. 

The inversion method used here is based on the regularized least-squares method- 
ology (RLS) following the prescription described in Eff-Darwich and Perez Hernandez 
(1997). In this implementation the regularization function is weighted differently for 
each model grid point (see also Eff-Darwich et al, 2007). In summary, Equation (2) 
is transformed into a matrix relation: 

D = Ax + e (3) 

where D is the data vector, with elements Avnim. and dimension M, x is the solution 
vector to be determined at A'' model grid points, A is the matrix with the kernels, 
of dimension M X N, and e is the vector containing the error bars. 
The RLS solution for the vector x is given by: 

xcst = {A'^ A + -rH)-^A'^D (4) 

where 7 is a scalar introduced to give a suitable weight to the constraint matrix H 
on the solution. Replacing D from Equation (3) we obtain 

xest = {A'^ A + -fHy^A'^ Ax = Rx (5) 

hence 

R= {A^A + jHy^A^A . (6) 

The matrix R is referred to as the resolution or sensitivity matrix. Ideally, R would 
be the identity matrix, which corresponds to perfect resolution. However, if we try to 
find an inverse with a resolution matrix R close to the identity matrix, the solution 
is generally dominated by noise magnification. The individual columns of R display 
how anomalies in the corresponding model are imaged by the combined effect of 
measurement and inversion. In this sense, each element Rij reveals how much of the 
anomaly in the j*^ inversion model grid point is transferred into the i*'' grid point. 
Consequently, the diagonal elements Ru indicate how much of the information is 
saved in the model estimate and may be interpreted as the resolvability or sensitivity 
of Xi. We defined the sensitivity \i of the grid point Xi to the inversion process as 
follow: 
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With this definition, a lower value of Ai means a lower sensitivity of Xi to the 
inversion of the solar rotation. We define a smoothing vector W with elements Wi = 
X~ that is introduced in Equation (4) to complement the smoothing parameter 7, 
namely 



Xest = (A^A + jWHy^A^D . (8) 

Such substitution allows us to apply difi^erent regularizations to difi^erent model 
grid points Xi, whose sensitivities depend on the data set that is used in the inver- 
sions. A set of results can be calculated for different values of 7, the optimal solution 
being the one with the best trade-off between error propagation and the quadratic 
difference = \Axcst—D\'^ as discussed in Eff-Darwich and Perez Hernandez (1997). 



Table 2. Description of the artificial data sets used to study 
the sensitivity of low-degree high-order p modes. 



Data set 


e=l (mHz) 


£ = 2,3 (niHz) 


e>3 (niHz) 


1 


1 < 1/ < 2.3 


1 < 1/ < 2.3 


1 < V < 3.9 


2 


1 < 1/ < 2.5 


1 < 1/ < 3.4 


1 < 1^ < 3.9 


3 


1 < < 3.4 


1 < < 3.4 


1 < V < 3.9 



We have performed a theoretical study to determine the effect of adding the low- 
degree high-order p modes in the inversions. To do so, we have computed, using 
Equation (2), the splittings corresponding to an artificial rotation profile {fl{r,9)). 
This artificial profile has a differential rotation in the convection zone (that mimics 
the real one), a rigid rotation from 0.7 down to 0.2 Rq equal to Qrz = 433 nHz and 
a step-like profile in the core having a rate of 350 nHz in the 0.1-0.2 Rq region and 
a rate three times larger than the rest of the radiative zone below 0.1 Rq. Although 
this profile with steep changes and a small drop followed by an increase is unlikely to 
be realistic, it enables us to characterize the quality of the inversions as the inversion 
code has difficulty reproducing these steep gradients. 

Three artificial data sets have been built (see Table 2). The first, Set 1, is the 
reference and contains low-degree p modes below 2.3 mHz while p modes for ^ > 3 
extend up to 3.9 mHz. The second set, Set 2, adds the 1 = 2 and 3 modes up to 
3.4 mHz. As we have shown in the previous section, these modes could have some 
sensitivity to the rotation in the core (see Figure 4a). Finally, the third set. Set 3, 
also adds the £ = I modes up to 3.4 mHz that seem to be at the sensitivity limit 
(see again Figure 4a). 

Figure 5 shows the resulting inversions for these three artificial data sets. We shall 
qualify these inversions as ideal because all of the splittings were given the same error 
bars and are noiseless. The inferred rotation rate represents the best result that we 
can obtain for each set of modes. 
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Figure 5. Equatorial rotation rate obtained from ideal inversions (all the modes with the same 
error bars) of an artificial step profile (continuous line) using the three data sets described in 
Table 2. The inversion error bars are plotted only for the profiles resulting from inverting Sets 
1 and 2. Shaded area corresponds to the region wrhere the solution of the inversion in not 
reliable as the result of the poor quality of the resolution kernels due to the lack of sensitivity 
of the acoustic modes to this region of the Sun. 



As expected, the differences between the profiles appear below 0.3 Rq but they are 
only significant below « 0.15 Rq. All along the solar interior, the differences between 
Set 2 and 3 are very small and are within the inversion uncertainties. Comparing 
these two data sets with Set 1, we obtain an improvement up to « 30% in the deepest 
region when the sets with more modes are used. Unfortunately, in all the data sets, 
the recovered profiles are not accurate enough. This is easily understood when looking 
at the corresponding resolution kernels: they remain, indeed, very broad. Figure 6 
shows some resolution kernels for the deepest model grid points and for Sets 1 and 
3. There are no resolution kernels centered below ~ O-lSii© for Set 1. Comparing 
the resolution kernels of Set 1 and 3, we see an improvement in the second case. 
Indeed the resolution kernels are narrower (better resolution in the solution of the 
inversion) and they are slightly shifted to deeper layers (by about 3% of the solar 
radius). 

We have shown the improvements of adding low-degree high-order p modes in 
the ideal inversions. However, the observed splittings have non-uniform error bars. 
Therefore, we inverted the same sets of artificial splittings but using the actual 
measured error bars. To be more realistic, Gaussian noise was also added to the 
artificial splittings, commensurate to their corresponding uncertainty. The resulting 
rotation profiles are plotted in Figure 7. As for the ideal case. Sets 2 and 3 give the 
same qualitative results and both are slightly better in the core than Set 1, even 
though the differences are within the inversion error bars. Nevertheless, the solution 
is improved in all of the radiative region where the ripples are reduced and the 
solution is smoother when more modes are added. Therefore, even with the realistic 
uncertainties, the recovered rotation rate is more accurate when Sets 2 or 3 are tised. 
Unfortunately, the contribution of the £ = I modes in the 2.5 to 3.4 mHz range 
remains negligible with the present magnitude of their uncertainties. 
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Figure 6. Resolution kernels of the inner model grid for Set 1 (top) and 3 (bottom) in the 
ideal case (all the modes have the same error bars). 
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Figure 7. Equatorial rotation rate obtained from the inversions of the artificial step profile 
(continuous line) using the three data sets described in Table 2 but with the real error bars 
as described in Section 2. For clarity, the inversion error bars are plotted only in the inversion 
of the Set 2. Shaded area corresponds to the region where the solution of the inversion in not 
reliable as explained in Figure 5. 



Including realistic uncertainties results in significant changes in the resolution 
kernels (see Figure 8) . Indeed, the coefficients of the linear combination given to each 
splitting (inversely proportional to the error bars) changes the linear combination of 
the modes used to compute those resolution kernels. Therefore the inner resolution 
kernels shift outwards: the maximum of the deepest one is now at around 0.16 Rq 
for all sets and, of course, are broader than the ones of the ideal case. The differences 
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between the resolution kernels deduced from Sets 1 and 3 are smaller, explaining 
why the solutions in both cases are so similar. 
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Figure 8. Resolution kernels of the inner model grid points for Set 1 (top) and 3 (bottom) 
but using splittings with realistic uncertainties. 

After having studied the improvements in the rotation rate that we can expect by 
adding the low-degree high-order p modes, we computed the same set of inversions 
using the real data. Figure 9 shows the inferred rotation rate for the real data using 
the Sets 1, 2, and 3. The results are similar to what we found for the artificial 
cases with realistic uncertainties. Indeed, the rotation rate below wO.25 Rq changes 
shghtly when comparing Set 2 and 3 (that cannot be distinguished within the error 
bars) to Set 1. The rotation profiles show a small increase inside the core when the 
new modes are added (Sets 2 and 3). The difference is barely significant down to 
0.16 Rq which is the reliability limit of the inversion. As we already knew, we need 
gravity modes to get more information about the dynamics of the deepest layers of 
the core (Mathur et al, 2007). In the rest of the radiative zone (above 0.25 Rq), the 
profiles inferred using the three sets are the same. 

We have also tested the effect of increasing the amplitude ratio of the visible m 
components of the 1 = 2 and 3 modes by 20 and 10% respectively. The resulting 
splittings are systematically increased by w 10 nHz for the £ = 2 modes and by 3 
to 4 nHz for the ^ = 3 modes. The inversion computed with these new splittings 
gives the same qualitative results as those shown in Figure 9: the inferred profiles 
are improved when more modes are used (Set 2 and 3). For these two profiles, the 
rotation rate obtained at 0.16 Rq has increased by some 7%, which is still within 
the inversion error bars. 



5. Conclusion 

We have analyzed a 4 182 day-long time series of GOLF velocity observation and 
have derived the frequencies and the splittings of the acoustic modes up to 3.4 mHz. 
We used a 45 /iHz window to fit the = 0, 2 and 1, 3 pairs to limit any bias in the 
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Figure 9. Equatorial rotation rate of real data using only low-degree low-order p modes up 
to 2.3 and 3.4 mHz and the same sot of high-degree modes. For clarity, the inversion error bars 
are shown only for the result of inverting Set 1. Shaded area corresponds to the region where 
the solution of the inversion in not reliable as explained in Figure 5. 

fitted splittings. By fitting shorter time series simultaneously and including a linear 
dependence of the splitting with the solar activity, we have checked that for most of 
the modes the variation of this parameter with the activity is within the error bars 
estimated using the whole time series. 

We then carried out a study of the sensitivity of the rotation of the core (below 
0.25 Rq) to the £ < 3 modes up to 3.4 mHz. We show that the contribution of the 
rotation rate in the core to the splittings of the £ = 2 and 3 modes is several times 
the size of their current error bars. By contrast the i = I modes in the 2.5 to 3.4 
mHz range have a contribution from the core of the same order of magnitude as 
their error bars. Therefore, they cannot provide valuable information on the core's 
rotation. 

By using artificial inversions of a step-like profile we have shown that, ideally, the 
new modes would improve the inferred rotation rate in the core but with a small 
accuracy on the profile itself due to a lack of resolution in the inner regions of the 
Sun - as indicated by the corresponding resolution kernels. The effect of including 
the £ = 1 modes between 2.5 and 3.4 mHz is negligible with the current size of their 
error bars. Using the real error bars degrades the inferred solution and the differences 
in the core profile are within the inversion uncertainties. However, the profile in the 
rest of the radiative region is smoother when more modes are used, demonstrating 
that the inclusion of those modes improves the inversion results. 

Finally, the inversion of the real data, including the larger set of low-degree low- 
order p modes, shows no significant differences between Set 2 and 3 as for the artificial 
data. There is a small increase of the rotation rate below 0.25 Rq compared to the 
reference Set 1 but it is within the inversion error bars. The resolution kernels are 
improved when the low-degree high-order p modes are included and the resolution 
kernels at the deepest locations are shifted by 3% of the solar radius towards the 
center. This result suggests that with the present error bars it is better to use such 
modes in the inversions (at least for the high order £ = 2 and 3 modes). 
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To progress in the inference of the rotation rate in the solar core we would need 
to measure low-degree modes: acoustic, mixed, and gravity modes that will bring a 
direct signature of the dynamics of the still hidden core of the Sun. 
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